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ABSTRACT 


As  the  costs  for  maintaining  computer  communication  networks  are  rapid!}’  rising, 
it  is  particularly  important  to  design  the  network  efficiently.  The  objective  of  this 
thesis  is  to  model  the  minimum  cost  design  of  digital  communication  networks  and 
propose  a  heuristical  solution  approach  to  the  formulated  model. 

The  minimum  cost  design  has  been  modeled  as  a  zero-one  integer  programming 
problem.  The  Lagrangian  relaxation  method  and  subgradient  optimization  procedure 
have  been  used  to  find  reasonably  good  feasible  solutions. 

Although  the  reliability  requirement  for  computer  communication  networks  is  as 
important  as  the  cost  factor,  only  the  cost  factor  is  considered  in  the  context  of  this 
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I.  INTRODUCTION 


Computer  communication  networks  are  becomming  a  vital  part  of  today’s  indus¬ 
trial,  govermental,  financial,  and  service  organizations.  The  efficiency  and  quality  of 
their  internail  operations  and  external  services  depends  upon  how  effective  their  com¬ 
puter  communication  networks  are.  Design  of  computer  communication  networks 
requires  the  consideration  of  a  number  of  factors  such  as  performance  constraints, 
reliabiiit)’  constraints,  and  costs.  As  the  costs  for  maintaining  computer  communi¬ 
cation  networks  are  rapidly  rising,  it  is  particularly  important  to  design  the  network 
efficienth’;  i.e.,  the  emphasis  on  keeping  the  costs  low  in  the  design  of  networks  is 
becoming  evident. 

Due  to  the  inherent  complexity  of  network  design  problems,  the  comprehensive 
optimal  design  of  efficient  networks  is  usually  infeasible  except  for  very  simple  net¬ 
works.  The  objective  of  this  thesis  is  to  model  the  minimum  cost  design  of  digital 
communication  networks  and  propose  a  heuristical  solution  approach  to  the  formu¬ 
lated  model. 

It  is  assumed  that  node  locations  are  fixed  by  the  required  placement  of  terminals 
and  network  facilities.  Therefore  network  design  is  predominately  a  question  of  where 
to  put  links  and  what  their  data  carrying  capacities  should  be.  For  a  given  set  of  n 
nodes,  there  are  n(n-l)/2  possible  links.  It  is  seldom  necessary  for  them  to  be  fully 
connected.  Consequently,  there  is  an  inherent  design  tradeoff  between  equipment 
simplicity  and  low  cost  with  fewer  links.  Although  the  reliability  requirement  for 
computer  communication  networks  is  as  important  as  the  cost  factor,  only  the  cost 
factor  is  considered  in  the  context  of  this  thesis. 

The  minimum  cost  design  has  been  modeled  as  a  zero-one  integer  programming 
problem.  The  Lagrangian  relaxation  method  and  subgradient  optimization  procedure 
have  been  used  to  find  reasonably  good  feasible  solutions. 

The  organization  of  this  thesis  is  as  follows:  Chapter  II  is  an  overview  of  digital 
transmission  facilities  and  describes  Tl,  T3  circuits  and  digital  cross-connect  systems. 
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Chapter  III  formulates  a  model  of  the  minimum  cost  design  problem  and  introduces 
a  Lagrangian  relaxation  method  along  with  a  subgrcidient  optimization  procedure. 
How  to  implement  this  technique  in  network  design  is  also  described.  Chapter  FV 
presents  the  results  of  computational  tests.  Finally,  conclusion  and  recommendations 
for  further  research  are  offered  in  Chapter  V. 
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II.  DIGITAL  TRANSMISSION  FACILITIES 


A.  Overview  of  Digital  Transmission 

Digital  transmission  teciiniques  were  introduced  in  the  Bell  System  in  the  early 
1960s  for  efEcient  transport  of  voice  signals.  At  that  time,  data  for  computer  com¬ 
munications  accounted  for  a  small  percentage  of  network  traffic.  Therefore,  the  basic 
digital  transmission  structure  was  centered  around  human  voice  communication. 

Although  a  plot  of  amplitude  versus  time  for  a  speech  sample  shows  significant 
high-frequenc}’  components,  a  successful  telephone  conversation  requires  less  than 
4KHz  of  audio  bandwidth.  After  being  processed  by  a  low-pass  filter,  a  speech  sample 
loses  its  high-frequency  components  and  is  ready  to  be  digitized. 

Since  the  maximum  analog  frequency  to  be  reproduced  is  4KHz,  an  8Kbit/s  sam¬ 
pling  rate  was  adopted  according  to  the  Nyquist  theorem.  By  sampling  the  voice 
signal  every  125  microseconds,  its  essential  information  is  extracted  from  the  ana¬ 
log  format  and  is  ready  for  digital  encoding.  This  pulse- amplitude  modulated(PAl.'I) 
signal  contains  all  the  information  in  the  original  signal  up  to  approximate!}'  4KHz. 

The  modulation  scheme  that  w'as  chosen  for  the  early  digital  transmission  stan¬ 
dards  is  the  eaisiest  to  implement.  The  pulse  code  modulation(PCM)  is  currently 
used  worldwide  for  a  standard  voice  communication.  Figure  1  shows  how  the  voice 
signal  is  encoded  to  digital  signal.  The  pulse  amplitude  modulation(PAM)  signal  is 
quantized  by  mapping  into  discrete  amplitude  levels,  each  with  i  unique  binary  code. 
Naturally,  additional  discrete  mappings  reduce  the  quantization  error  and  improve 
the  fidelity  of  the  coded  signal.  The  analog  voice  signal  is  sampled  8000  times  per 
second  with  each  sample  measuring  the  amplitude  of  the  signal,  and  sampled  ampli¬ 
tude  is  encoded  into  an  8-bit  digital  code.  By  creating  8-bit  codes  8000  times  per 
second,  the  rate  of  a  digitized  voice  signd  is  8  multiplied  by  8000,  or  64Kbit/s.  This 
forms  the  b^Lsic  64Kbit/s  channel  that  is  the  foundation  of  the  digital  network.  Such 
a  digital  channel  is  often  referred  to  as  a  DS-0(digital  signal  level  0). 
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Figure  1:  Pulse  Code  Modulation 


The  original  deploymant  of  digital  transmission  wais  driven  by  a  desire  to  conserve 
copper  pairs  outside  telephone  offices.  Each  analog  voice  signal  in  the  existing  tele¬ 
phone  network  consumed  a  physical  pair  of  copper  wires  from  a  subscriber  location 
to  a  telephone  central  office.  Once  the  basic  block  was  digital,  it  became  feasible  to 
multiplex  the  digital  signals  together  into  a  higher-order  digital  signal.  The  time- 
division  multiplexing  led  to  defining  the  next  level  of  the  digital  hierarchy  as  being 
equivalent  to  24  DS-0  signals.  This  level  is  referred  to  as  DS-1  or  Tl.  [Ref. 13] 

The  multiplexing  function  that  created  a  DS-1  signal  was  originally  handled  by  3 
network  element  known  as  a  digital  channel  bank  or  D-bank.  Numerous  devices,  of 
course,  offer  Tl  interfaces  at  these  days. 

B.  Tl  System 
1.  What  is  Tl? 

The  Tl  system  has  become  one  of  the  most  widely  used  high-capacity  sys¬ 
tems  for  transmission  of  voice  and  data.  Originally  conceived  as  voice  transmission 
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technology  in  the  early  1960s,  it  has  evolved  into  a  cost-effective  and  highly  flexible 
means  of  transmitting  both  voice  and  data. 

Tl  is  based  on  time-division  multiplexing  (TDM)  24  users  onto  one  phys¬ 
ical  circuit.  It  usually  involves  a  digital  service  unit(DSU)  and  a  channel  service 
unit(CSU)  to  accomplish  the  analog-to-digital  conversion  and  the  multiplexing.  At 
the  inception  of  Tl,  it  was  known  that  solid  copper  cable  was  capable  of  provid¬ 
ing  frequencies  about  lOOKHz.  Consequently,  the  earlier  frequency  division  multi- 
plexed(FDM)  telephone  system  multiplexed  24  users  onto  one  copper  wire.  Each 
user  was  assigned  a  different  4-KHz  frequency  band.  The  composite  of  24-voice  chan¬ 
nels  totalled  96KHz,  which  was  within  the  capacity  of  twisted  pair  wires.  [Ref .11] 

The  Tl  system  is  designed  around  a  1.544Mbit/s  rate,  which  was  about  the 
highest  rate  that  could  be  supported  across  twisted  pair  for  a  distance  of  approxi¬ 
mately  one  mile. 

Tl  has  several  advantages.  The  economic  incentive  to  use  Tl  can  be  enormous 
despite  the  apparent  high  initial  cost  of  equipment  and  the  continuing  costs  for  leased 
lines.  An  investment  in  a  new  network  based  on  Tl,  however,  can  break  even  within 
one  year.  The  end  user  organization  operating  a  private  Tl  network  enjoys  the 
diagnostic  power  to  identify  and  isolate  faults  and  the  ability  to  set  up  and  change 
connections  between  points  on  the  network  at  will.  Another  advantage  is  that  Tl 
technology  simplifies,  and  thereby  improves  corporate  communications.  Combining 
multiple  voice  and  data  channels  over  a  single  high-speed  digital  circuit  eliminates 
many  separate  networks  and  many  forms  of  network  management.  Therefore,  Tl 
makes  network  control  more  practical  and  improves  reliability.  With  the  availability 
of  large  amounts  of  bandwidth,  under  immediate  control,  the  network  operator  is  able 
to  offer  the  corporation  new  services  that  were  not  practical  in  the  pjist.  Additional 
circuits  are  readily  available  through  reconfiguration.  [Ref.W] 

2.  Tl  Frame  Format 

The  majority  of  Tl  offerings  digitize  the  voice  signal  through  PCM  or  adaptive 
differential  pulse  code  modulation(ADPCM).  Whatever  the  encoding  technique  used. 
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Figure  2:  The  Tl  Based  Frame 

once  the  analog  images  are  translated  to  digital  bit  streams,  man\'  Tl  systems  then 
time-division  multiplex  voice  and  data  together  in  24  or  44  user  slots  within  each 
frame. 

This  section  introduces  the  concept  that  the  carriers  transmit  24  voice  or  data 
channels  together  with  TDM  frames.  The  Tl  carrier  system  provides  multiplexing 
by  sampling  the  24  channels  at  a  combined  rate  of  192.000  times  per  second.  Figure 
2  shows  how  the  24  channels  are  multiplexed  into  a  frame.  The  frame  contains  one 
sample  from  each  channel,  plus  an  additional  bit  for  frame  synchronization.  Thus, 
the  complete  frame  contains  193  bits  (8  bits  per  channel  x  24  channels  -I-  1  sync, 
bit  =  193bits).  In  addition,  each  of  the  24  channels  requires  a  64Kbit/s  rate  (8000 
samples/sx  8  bit/sample  =  64Kbit/s). 

The  baisic  Tl  frame  consist  of  24  eight-bit  slots  and  a  framing  bit.  The  details 
of  Tl  frame  format  are  characterized  by  the  following  concepts: 


•  Alternate  Mark  Inversion  (AMI) 

•  Is  density  requirements 

•  Bit  robbing 

•  Binary  7  zero  substitution  (B7ZS) 
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Figure  3:  Format  errors  Due  to  Bipolar  Code  Violation 
e  Framing  bit 


Alternate  Mark  Inversion 

The  AMI  concept  is  shown  in  Figure  3(a).  The  user  may  recognize  this  code  as 
bipolar  coding.  The  Tl  system  requires  that  a  1  pulse  must  be  sent  as  an  opposite 
polarity  from  the  preceding  1  pulse,  regardless  of  the  number  of  Os  in  between  the 
two  Is.  This  bipolar  code  performs  well  becavoe  it  hais  no  direct  current  component 
and  therefore  can  be  coupled  with  transformers.  The  scheme  also  makes  efficient  use 
of  bandwidth. 

Is  Density  Requirements 
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T1  provides  no  separate  clocking  signal.  The  timing  or  clocking  information  is 
embedded  in  the  data  stream.  At  the  receiving  end,  the  clock  is  recovered  from  the 
data  stream  by  the  detection  of  1  pulse.  If  the  Tl  data  stream  has  insufficient  1 
pulses  embedded,  the  receiver  can  no  longer  produce  reliable  timing  out. 

To  overcome  this  problem,  a  certain  number  of  Is  must  be  present  to  ensure  proper 
timing.  This  concept  is  called  Is  density.  The  Tl  facilities  require  that  no  more  than 
15  consecutive  Os  shall  be  present  in  the  frame.  An  additional  requirement  is  that 
there  must  be  at  least  three  Is  in  every  24  bits. 

Due  to  noise,  mechanical  failures,  etc.,  the  bits  may  become  distorted,  which  can 
cause  a  violation  of  the  AMI  rule.  Bipolar  violations  or  excessive  errors  are  known  as 
format  errors  because  the  errors  are  not  in  conformance  with  the  required  Tl  format. 
A  bit  distortion  may  not  cause  a  format  error.  As  Figure  3(b)  through  3(e)  show,  the 
nature  of  the  error  determines  if  the  bipolar  violation  is  detected.  The  bipolar  signal 
is  altered  with  errors  of  omission  or  errors  of  commission. 

Bit  Robbing 

The  least  significant  bit  of  every  time  slot  in  the  ’signalling  frames’  is  overwritten 
to  encode  the  control  signals.  This  concept  is  called  bit  robbing. 

For  the  transmission  of  data,  the  8th  bit  is  unavailable.  Consequently,  the  majority 
of  Tl  and  related  systems  use  a  56Kbit/s  transmission  rate  instead  of  the  64Kbit/s 
rate  actually  available. 

Binary  7  Zero  Substitution 

To  establish  Is  density  in  the  Tl  signal,  a  technique  called  B7  zero  code  suppres¬ 
sion  is  used.  As  previously  discussed,  a  1  is  substituted  in  a  Tl  frame  to  prevent  an 
occurrence  of  more  than  15  consecutive  Os.  It  is  possible  that  this  substitution  can 
occur  with  any  bit  in  the  frame.  A  voice  channel  can  tolerate  the  loss  of  the  7th  bit 
to  B7  substitution.  However,  a  data  channel  cannot  operate  with  this  arrangement. 

It  is  not  necessary  to  rob  bits  in  a  data  channel  since  there  is  no  telephone  signal. 
However,  it  is  quite  possible  to  have  all  Os  in  a  data  channel.  Therefore,  the  data 
could  be  corrupted  by  the  B7  zero  code  suppression  technique.  The  insertion  does 


8 


T-1  Facility  Termination 


f 

I 

I 

I 


DSX-i 

Inwriac* 


r 

Tl  Mux 

L_ 

11 

Vanous  ImertacM 

And  Connectors  V35.  RS-422 

Internal  to  Mux 


V 

+ 

D 


From  .  I  I  ^  - 1 !.  Cuttotnar  Suppbad 

Present 

Oemark 


Figure  4:  T1  Facility  Termination 


not  affect  the  voice  channel,  but  it  certainly  affects  a  data  channel.  To  prevent  this 
problem,  a  data  channel  generally  only  uses  7  bits  and  any  1  substitution  can  be 
placed  into  the  8th  bit.  This  prevents  the  B7  zero  code  suppression  from  corrupting 
the  data  channel. 

Framing  Bit 

To  decode  the  incoming  data  stream,  a  receiver  must  be  able  to  associate  each 
sample  with  the  proper  channel.  At  a  minimum,  the  begining  and  ending  of  the  frame 
must  be  recognized.  The  function  of  the  framing  bit  provides  this  requirement.  This 
bit  is  located  in  the  193rd  bit  of  each  frame.  It  is  not  part  of  user’s  information, 
but  added  by  the  system  for  framing.  The  use  of  framing  bit  varies  significantly 
depending  on  the  type  of  Tl  and  of  the  technology.  [Ref  .17] 


3.  Tl  Networks 

(A)  Tl  Interface 

Figure  4  shows  Tl  facility  termination.  The  Tl  facility  termination  enters 
the  users  premises  as  two  wire  pairs  for  data.  There  may  be  additional  wires  for 
loopback  control  and  testing  of  unframed  lines.  The  channel  service  unit  (CSU) 
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T-1  Trunks  Trunks 


Figure  5:  Digital  Cross-Connect  System 

provides  loopbacks  to  the  network  and  is  the  last  regeneration  point  for  incoming 
signals.  DSX-1  as  an  interface  is  well  defined  in  terms  of  electrical  and  mechanical 
features.  A  digital  service  unit  (DSU)  matches  various  terminal  interfaces  to  the 
DSX-1  standard  and  usually  is  part  of  the  user’s  terminal  equipment:  multiplexer, 
bandwidth  manager,  PBX,  computer,  etc.[/?e/.ll] 

(B)  Digital  Cross-Connect  Systems 

Digital  cross-connect  systems  (DCS)  are  electronic  switches  able  to  take 
bit  streams  out  of  the  Tl-carrier  and  route  them  to  a  given  output  line.  DCS  allow 
remote  access  to  and  cross-connection  of  large  bandwidth  circuits.  Figure  5  shows 
DCS.  DCS  can  be  divided  into  three  primary  subsystems. 

At  the  core  of  the  DCS  is  the  Matrix  Subsystem  where  the  actual  chan¬ 
nel  level  (DS-0)  cross-connections  take  place.  DCS  systems  are  controlled  by  the 
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Data  Bypass 

Node  B 


Figure  6:  Data  Bypass 

Administrative  Subsystem,  where  high-level  processors  manage  the  oveiall  system 
functions.  Other  common  elements  such  as  system  synchronization  and  user  inter¬ 
face/control  port  access  reside  within  this  subsystem.  Finally,  there  is  the  Interface 
Subsystem,  which  terminates  DS-1  or  higher  facilities.  Interface  subsystems  are  orga¬ 
nized  into  growth  increments,  or  units,  which  accommodate  anywhere  from  28  to  32 
slot  positions  each.  Aside  from  facility  terminations,  these  slot  positions  can  be  used 
to  accommodate  features  such  as  sub-rate  date  multiplexing  and  bridging  of  digital 
services. 

Figure  6  shows  data  bypass  based  on  DCS.  Data  bypass  resembles  the 
drop  and  insert  function  of  the  voice  channel  bank.  What  distinguishes  the  two  is 
that  data  bypass  offers  variable  bandwidth  allocation  and  additional  flexibility  in 
routing  connection  among  many  more  Tl  links. [ile/.  11] 

(C)  Tl  Networking 

Lower  cost  for  digital  lines  and  for  customer  premises  equipment  have 
made  all-digital  networking  not  only  possible,  but  attractive.  Every  location  is  on 
the  digital  backbone,  through  small  sites  may  be  served  by  56Kbit/s  or  other  data 
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Figure  7:  Tl  Networking 

link  speed  less  than  a  full  Tl. 

The  Tl  networking  are  shown  in  Figure  7.  There  are  two  main  benefits 
in  Tl  networking.  Trends  in  tariffed  costs  are  up  for  analog  facilities  and  down  for 
digital  facilities.  As  time  goes  on,  it  takes  fewer  analog  lines  to  justify  a  digital  line. 
Currently,  as  few  as  3  or  4  leased  voice  grade  lines  cost  as  much  as  a  56kbit/s  digital 
service.  As  few  as  two  or  three  56K  lines  cost  as  much  as  a  Tl.  The  result  is  that  we 
can  cost-justify  digital  service  in  more  and  more  locations.  Although  there  exist  var¬ 
ious  types  of  nodes  with  a  common  network  management  system,  everything  can  be 
brought  under  one  network  management  system  because  Tl  networking  is  all-digital. 
[Ref.n] 
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C.  T3  SYSTEM 

1.  What  is  T3? 

Many  corporate  managers  are  administering  T1  networks  that  are  growing 
rapidly  with  new  users  asking  for  more  amd  more  bandwidth.  Increasingly,  data 
communications  managers  are  eyeing  T3  network  as  a  solution  to  the  problems  raised 
by  this  expansion. 

T3  refers  to  a  digital  transmission  facility  running  at  DS-3  data  rate  of 
44.736Mbit/s.  It  is  equivalent  of  28  Tls,  and  is  more  economical  than  a  much  less 
number  of  Tls.  For  a  distance  of  less  than  50  miles,  only  four  Tls  are  required  to 
break  even. 

T3  circuits  usually  use  microwave  and  optical  fiber  as  transmission  media. 
Optical  fiber  has  a  much  higher  transmission  capacity  than  digital  microwave.  Dig¬ 
ital  microwave  has  some  problems.  First,  the  electromagnetic  spectrum  is  strictly 
licensed  by  the  Federal  Communications  Commission.  Second,  strict  spectrum  al¬ 
locations  imply  that  digital  microwave  cannot  be  readily  upgraded  in  bandwidth. 
While  a  45Mbit/s  optical  fiber  can  be  readily  converted  to  565Mbit/s,  a  45Mbit/s 
digital  microwave  setup  is  probably  destined  to  remain  at  45Mbit/s.  Finally,  digital 
microwave  devices  are  inherently  limited  to  line-of-sight  distance. 

Many  T3  networks  mix  optical  fiber  and  digital  microwave.  The  optical  fiber 
can  be  used  in  high-density  areas,  with  microwave  hops  to  lower-density  remote  lo¬ 
cations.  Alternatively,  the  main  network  can  be  placed  entirely  on  fiber,  with  T3 
microwave  used  for  emergency  restoration  in  the  event  of  a  cable  cut.  Dual-media 
transmission  is  especially  attractive  to  service-critical  networks  such  as  banks,  airlines, 
and  other  on-line  transaction-processing  companies. 

T3  was  originally  designed  to  be  used  as  interexchange  trunk  within  the  public 
telephone  network. [i?e/.  13] 

2.  T3  Format 

By  the  end  of  the  1970s,  a  way  around  the  limitations  of  the  copper  pair 
was  on  the  horizon.  Transmission  technologies  baised  on  the  nearly  limitless  band- 


13 


width  of  optical  fiber  were  emerging  from  the  laboratory,  and  a  new  layer  of  the 
digital  hierarchy  seemed  appropriate.  Based  on  technology  available  at  the  time, 
an  asynchronous  DS-3  rate  was  defined  as  the  combination  of  seven  DS-2  signals  of 
6.176Mbit/s.  Framing  and  stuffing  bits,  as  well  as  rudimentary  error  checking  and 
internal  communications,  were  added  to  the  DS-2  tributaries  to  generate  an  aggregate 
bit  rate  of  44.736Mbit/s. 

Notice  that  the  information  content,  or  payload,  of  a  DS-3  is  equal  to  672 
times  64Kbit/s  =  43.008Mbit/s.  The  additional  1.728Mbit/s  of  overhead  represents 
the  sum  of  the  DS-1  framing  bit,  the  DS-2  framing  and  stuffing  bits,  and  the  variety  of 
DS-3  overhead  bits.  In  terms  of  network  efficieny,  the  DS-3  format  devotes  96  percent 
of  the  transmission  bandwidth  to  payload,  with  approximately  4  percent  overhead. 
[i2e/.13] 

Although  the  DS-3  was  designed  to  be  created  from  multiple  DS-2  signals, 
there  was  an  obvious  inefficiency  in  using  two  separate  devices  for  the  DS-l-to-DS-2 
and  DS-2-to-DS-3  multiplexing  functions. 

There  are  few  choices  when  it  comes  to  T3  formats.  The  long  established 
standard  in  the  public  network  is  now  called  M13.  But  for  technical  reasons  there 
are  new  proposals  for  formats  called  Syntran  and  SONET. 

(A)  Ml3  Format 

The  name  M13  comes  from  the  fact  that  this  multiplexer  connects  the 
DS-1  to  the  DS-3  level.  One  DS-3  contains  28  DS-ls.  They  are  combined  in  two  steps. 
The  M13  process  first  joins  four  DS-1  lines  into  DS-2  bit  stream  at  6.312Mbit/s.  This 
step  is  more  than  a  straight  time  division  multiplexer  function.  M12  bit-interleaves 
the  four  inputs,  but  it  also  adds  136Kbit/s  of  overhead  and  justification  or  bit  stuffing. 

The  second  step  within  M13  combines  seven  DS-2  streams  into  one  DS-3. 
Like  M12,  the  M23  process  is  bit-interleaved,  with  justification.  That  is,  the  DS-3 
clock  need  not  be  synchronized  with  any  DS-2  or  DS-1  line  speed.  Ml 3  overruns  the 
DS-2  clock,  and  bit  stuffs  or  adds  extra  bits  to  fill  out  the  DS-3  rate  of  44.736Mbit/s. 

Whatever  the  original  reasons  for  the  two-step  multiplexing,  it  creats  a 
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problem  for  us  today.  There  is  no  reliable  way  to  tell  where  a  given  DS-0  is  located 
within  the  DS-3  bit  stream.  This  makes  it  very  difficult  to  extract  just  one  DS-0 
because  it  will  have  no  definite  location  within  the  DS-3  frame. 

To  get  at  an  individual  channel  you  must  typically  reverse  the  two  step 
and  demux  the  DS-3  down  to  individual  DS-ls.  The  need  for  a  DS-3  digital  access 
and  cross-connect  system(DACS)  is  clear.  To  get  simple  switching  capability  among 
DS-3  lines,  a  better  format  is  required.  [Ref.ll] 

(B)  Syntran 

The  Syntran  has  been  proposed  by  Bellcore.  Syntran  describes  how  to 
multiplex  DS-0  and  DS-1  signals  into  a  DS-3,  in  a  fully  synchronous  format. 

The  resulting  bit  stream  has  the  same  rate  as  the  standard.  Syntran  offers 
two  important  features. 

•  Synchronous  format  allows  easy  DACS  switching  of  DS-0  and  DS-1  channels 
among  DS-3  ports. 

•  The  aggregate  bit  rate  is  compatible  with  the  large  number  of  existing  DS-3 
transmission  facilities  installed  in  North  America. 

Unfortunately  for  the  global  network,  DS-3  is  seldom  found  outside  of 
North  America.  And  while  M13  and  Syntran  are  defined  as  electrical  interfaces,  the 
new  installations  of  facilities  near  that  speed  are  mostly  optical  fiber.[ile/.13] 

(C)  SONET 

SONET(Synchronous  Optical  Network)  is  much  more  than  DS-3.  The 
goal  is  a  standard  for  interconnection  of  different  national  networks  at  speeds  from 
about  150Mbit/s  up  to  several  gigabits/s.  The  DS-3  level  is  only  one  portion  of  the 
standard,  intended  for  network  access  from  customer  premises. 

SONET  operates  at  multiples  of  T3  bandwidth.  Initially,  products  will  be 
offered  at  the  following  bit  rates  (OC  stands  for  Optical  Carrier):  OC-1,  51.84Mbit/s; 
OC-3,  155.52Mbit/s;  OC-12,  622.08Mbit/s;  OC-48,  2.49Gbit/s. 
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As  a  synchronous  standard,  SONET  is  well-suited  for  switching  tributary 
signals  within  a  higher  bandwidth  pipe.  Initially,  switching  will  be  limited  to  T1  and 
DS-0  signals,  but  other  service  offerings  will  be  defined  in  the  future. 

SONET  is  fully  backward- compatible  with  the  ANSI  synchronous  T3  for¬ 
mat;  it  is  partially  backward-compatible  with  order  asynchronous  T3  equipment  such 
as  Ml3s.  It  does  not  affect  users  of  digital  microwave  radio,  since  there  is  little  need 
for  a  synchronous  standard  higher  than  T3  for  these  applications.  [Ref. 13] 

3.  T3  Networks 

Implementing  a  T3  corporate  backbone  provides  users  with  a  measure  of  flex¬ 
ibility  and  control  over  the  network.  At  the  low-speed  end,  circuits  may  be  either 
56Kbit/s  or  64Kbit/s[ile/.13].  These  may  be  used  for  voice  or  data  services. 

The  capability  to  switch  the  variety  of  circuits  allows  the  user  to  put  up  and 
tadce  down  circuits  of  differing  rates  on  the  T3  backbone  as  requirements  change, 
without  having  to  work  through  an  external  network  provider.  This  can  take  the 
form  of  reducing  lead  time  for  a  circuit  order,  time-of-day  circuit  changes,  or  quick 
response  to  a  temporary  overload  condition  in  some  part  of  the  network. 

T3  networks  are  based  on  one  or  more  of  the  following  types  of  DOS’s. 

1.  Point  to  point  terminals:  M13  multiplexers  with  28  Tls  in  and  one  DS-3 
out. 

2.  Multiple  Ml 3  multiplexers  in  a  single  chassis:  something  like  a  48-port  D4 
channel  bank  that  followed  the  24-port  D3. 

3.  Networking  of  multiple  DS-3  lines  into  one  hub  with  the  ability  to  cross- 
connect  DS-ls  and  even  DS-Os  within  that  hub. 

Level  3  generally  is  possible  only  with  one  of  the  synchronous  formats,  not 
with  classic  M13.  Devices  with  this  performance  are  installed  at  operating  companies. 

Network  topologies  for  T3  backbones  will  evolve  from  point  to  point  up  to 
the  partial  mesh  or  interlocking  rings  familiar  from  Tl  nets.  There  will  be  many 
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fewer  sites  that  cam  cost-justify  multiple  DS-3s  than  there  wih  be  Tl  hubs.  The 
combination  likely  to  be  most  common  is  a  sparse  DS-3  backbone  surrounded  by  a 
Tl  network.  [J?e/.ll] 
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m.  MODEL  FORMULATION  AND  SOLUTION 


This  chapter  discusses  how  to  formulate  the  minimum  cost  design  of  T1/T3  net¬ 
works,  and  explains  the  heuristic  procedure  for  optimal  solutions  based  on  the  La- 
grangian  relaxation  method. 

A.  MATHEMATICAL  FORMULATION 

In  order  to  model  the  minimum  cost  design  of  a  network,  we  consider  several 
things.  The  amount  of  traffic  between  each  node  and  every  other  node  is  given  and 
we  define  the  set  of  circuit,  types  each  of  which  is  a  mix  of  Tl ,  T3  and  fractional  Tl .  If  a 
network  contains  n  nodes,  there  would  be  n(n-2)/2  communicating  source-destination 
pairs  that  select  routes  from  sets  containing  from  a  few  thousand  to  hundreds  of 
thousands  of  candidate  routes,  leading  to  very  large  problem  sizes.  Therefore,  to 
simplify  the  computation,  we  assume  that  the  maximum  length  of  the  path  for  each 
pair  is  less  than  a  specified  number  of  hops.  For  example,  the  network  designer 
may  specify  that  routes  between  two  nodes  should  peiss  through  no  more  than  two 
additional  nodes. 

We  introduce  the  following  notations  to  present  the  mathematical  formulation  of 
the  problem. 

P  :  The  index  set  of  the  communicating  source-destination  pairs 

Rp  ;  The  index  set  of  candidate  routes  for  the  source-destination  pair  p  E  P.  The 
set  of  candidate  routes  is  generated  by  a  route  generation  algorithm.  A  route 
is  characterized  by  the  ordered  set  of  links  from  source  to  destination  node. 

L  :  The  index  set  of  candidate  links  in  the  T1/T3  network. 

T  :  The  index  set  of  circuit  types. 

Xpr  :  A  decision  variable  which  is  one  if  the  source-destination  pair  p  takes  the  route 
r  and  zero  otherwise. 
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Zit  :  A  decision  variable  which  is  one  if  the  link  /  is  supported  by  the  circuit  type  t 
and  zero  otherwise. 

Qt  :  The  number  of  DS-0  channels  that  can  be  carried  over  a  circuit  type  t. 
ap  :  The  number  of  DS-0  channels  demanded  by  the  source-destination  pair  p. 

Cit  :  The  costs  for  leasing  a  circuit  type  t  on  link  1. 

Iri  :  An  indicator  function  which  is  one  if  link  I  is  used  in  route  r  and  zero  otherwise. 

Using  the  above  notations,  we  formulate  the  minimum  cost  design  problem  as 
follows; 


Problem  P: 
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The  constraint  (2)  ensures  that  each  source-destination  pair  selects  exactly  one 
route  for  transmission.  The  constraint  (3)  ensures  that  the  number  of  DS-0  channels 
on  each  link  does  not  exceed  its  capacity.  The  selection  of  at  most  one  circuit  type 
for  each  link  is  ensured  by  constraint  (4).  Problem  F  is  a  linear  zero-one  integer 
programming  problem  which  can  be  shown  to  be  NP-hard. 

Based  on  the  current  status  of  the  computational  ability,  the  idea  of  solving 
Problem  P  optimally  is  rejected.  The  idea  of  using  piecewise  linearization  of  the 
objective  function  and  solving  the  linear  programming  problem  is  also  rejected  for 
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the  reason  that  the  linear  relaxation  of  Problem  P  is  very  large  and  its  optimal  solu¬ 
tion  requires  a  significant  coniputational  effort  and  may  generate  a  fractional  solution. 

B.  LAGRANGIAN  RELAXATION 

of  the  most  computationally  useful  ideas  of  the  1970s  is  the  observation  that 
many  hard  problems  can  be  viewed  as  easy  problems  complicated  by  a  relatively  small 
set  of  side  constraints.  Dualizing  the  side  constraints  produces  a  Lagrangicin  problem 
that  is  easy  to  solve  and  whose  optimal  value  is  a  lower  bound  for  minimization 
problems  on  the  optimal  value  of  the  original  problem.  The  Lagrangian  problem 
can  thus  be  used  in  place  of  a  linear  programming  relcixation  to  provide  bounds  in  a 
branch  and  bound  algorithm.  The  Lagrangian  approach  offers  a  number  of  important 
advcintages  over  linear  relaxation. 

Lagrangian  methods  had  gained  considerable  currency  by  the  1974  when  Geoffrion 
coined  the  perfect  name  for  this  approach  -  Lagrangian  relaxation.  Since  then  the 
list  of  applications  of  Lagrangian  relaxation  has  grown  to  include  over  a  dozen  of  the 
most  infamous  combinatorial  optimization  problems.  For  most  of  these  problems, 
Lagrangian  relaxation  has  provided  the  best  existing  algorithm  and  has  enabled  the 
solution  of  problems  of  practical  size.  [Ref  .3] 

In  the  previous  section,  we  already  made  a  mathematical  formulation  for  minimum 
cost  design  problem.  We  know  that  there  are  two  natural  Lagrangian  relaxation  for 
Problem  P.  The  first  is  obtained  by  multipling  the  constraint  (3)  by  a  vector  of 
positive  Lagrange  multipliers,  Ui,  1  €  L,  and  adding  them  to  the  objective  function. 

Problem  D\ : 
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The  set  of  feasible  solutions  for  Problem  F  is  a  subset  of  the  set  of  feasible  solutions 
for  the  Lagrangian  relaxation  Di.  The  positivity  of  U  ensures  that  in  any  feasible 
solution  for  Problem  P,  the  expression 

I€L  ipePreRp  teT 

is  nonpositive,  and  thus,  the  value  of  the  objective  function  in  (6)  is  never  greater 
than  the  value  of  the  objective  function  in  Problem  P.  Thus,  whenever  Problem 
P  has  a  feasible  solution,  Zoi{U)  is  less  than  or  equal  to  Zp.  For  each  vector  of 
multipliers  [/,  Zd\{U)  is  a  lower  bound  for  Zp. 

The  second  relaxation  is  obtained  by  dualizing  constraint  (2)  with  nonpositive 
Lagrange  multipliers, Vp,  p  ^  P. 

Problem  D^: 

ZB2{v)  =  wiu^^CuZH+'£K  EAV-1  (7) 
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Comparing  above  two  Lagrangian  relaxations,  we  see  that  the  second  relaxation  is 
harder  to  solve  but  might  provide  better  bounds.  This  thesis  uses  the  first  Lagrangian 
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relaxation  for  implementation  and  computational  results. 

C.  SUBGRADIENT  OPTIMIZATION  PROCEDURE 

Several  methods  have  been  suggested  in  the  hterature  for  computing  a  vector  of 
optimal  Lagrange  multipliers.  These  include  various  versions  of  the  simplex  method 
implemented  using  column  generation  procedures,  dual  ascent  and  multiplier  ad¬ 
justment  procedures,  and  a  subgradient  optimization  procedure[ile/.3].  Subgradient 
optimization  methods  have  been  shown  to  be  very  effective  in  a  variety  of  combinato¬ 
rial  optimization  problems  such  as  the  traveling  salesman  problem,  topological  design 
of  computer  communication  systems,  and  lot-sizing  algorithms  for  complex  product 
structures[i?e/.l].  Therefore,  we  decided  to  employ  subgradient  optimization  tech¬ 
niques. 

The  subgradient  method  is  a  brazen  adaptation  of  the  gradient  method  in  which 
gradients  are  replaced  by  subgradients. 

Let  (Xpr{U),  Zit{U))  be  the  optimal  solution  to  the  Lagrangian  problem  D\  for  a 
fixed  vector  U.  A  subgradient  is  the  vector  with  coordinates 

-n(C')  =  E  E  -  E 

peP  r€ftp  ter 

Given  an  initial  value  t/°,  a  sequence  is  generated  by  the  iterative  formula 

=  Ut  +  hit 

where  tk  is  a  positive  scalar  stepsize.  It  is  known  that  17*  converges  to  the  optimal 
Lagrange  multiplierL^’  provided  that  tk  converges  to  0  and  h  diverges.  The  stepsize 
most  commonly  used  in  practice  is 

A»  [Zp  -  Zp,(l/^)] 

'*  “  II  -if  r 

where  is  a  scalar  satisfying  0  <  A*  <  2  and  Zp  is  an  upper  bound  on  Zxjj(L), 
frequently  obtained  by  applying  a  heuristic  to  Problem  P.  Often  the  sequence  A*  is 
determined  by  setting  A©  =  2  and  halving  A*  whenever  Zd\{U)  has  failed  to  increase 
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in  some  fixed  number  of  iterations.  This  rule  has  performed  well  empirically,  even 
though  it  is  not  guaranteed  to  satisfy  the  sufficient  condition  given  above  for  optimal 
convergence. 

The  steps  involved  in  the  subgradient  optimization  procedure  are  as  follows. 

1.  Initialization: 

(a)  Using  a  heuristic,  compute  an  overestimate  Zp  of  Zp  or  set  Zp  to  an 
arbitrarily  large  value. 

(b)  Select  an  initial  set  of  multipliers  and  set  K,  the  iteration  counter,  to 
0,  and  the  improvement  counter  to  0;  and  A  a  parameter  for  adjusting  the 
stepsize  to  A°  (an  arbitrary  positive  initial  value,  e.g.,  2). 

2.  Solving  the  Lagrangian  problem: 

(a)  Increment-  the  improvement  and  iteration  counters  by  1. 

(b)  Solve  the  Lagrangian  problem  using  as  the  Lagrange  multipliers.  Thus, 
obtain  Zdi{U‘‘)  and  Zf^. 

3.  Testing  and  updating  parameters: 

(a)  If  Zo\{U^)  is  greater  than  the  current  best  value  of  Zdi{U)  then  replace 
the  current  best  value  of  Zdi{U)  by  Zdi{U^)  and  set  U‘  to  U‘‘.  Also  reset 
the  improvement  counter  to  -1. 

(b)  If  and  Zj^  are  feasible  for  Problem  P,  compute  the  value  of  its  as¬ 
sociated  object! velfunction  for  Problem  P.  If  this  value  is  less  than  the 
current  value  of  Zp,  then  set  Zp  to  this  value. 

(c)  If  the  improvement  counter  has  reached  a  prespecified  upper  hmit,  then 

set  A  to  A/2,  to^*,  the  improvement  counter  to  0  and  go  to  2(a). 

(d)  If  the  iteration  counter  has  exceeded  a  prespecified  limit,  if  A  is  less  than 
a  prespecified  limit,  if  t*  is  less  than  a  prespecified  limit,  or  if  {Zp  — 
Zd\{^’)) I Zd\{U')  is  less  than  a  prespecified  error  tolerance,  then  stop. 
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4.  Updating  the  multipliers  : 

(a)  Compute  a  new  subgradient. 

=  E  E  -  E  QX  «  e  i 

p€P  reRp  teT 

(b)  Compute  the  new  stepsize. 

A[Zp-ZDi(y*)] 

‘‘  ■  li  7?  r 

(c)  Compute  the  new  multipliers. 

=Ut  +  tklt  V/€l 

(d)  set  k  to  k  +  1. 

5.  Go  to  2. 

D.  SOLVING  THE  LAGRANGIAN  PROBLEM 

The  computational  efficiency  of  the  subgradient  optimization  procedure  depends 
on  our  ability  to  solve  efficiently  the  Lagrangian  problem  which  is  generated  in  step 
2(b)  of  the  subgradient  optimization  procedure.  Fortunately,  for  a  fixed  set  of  multi¬ 
pliers,  the  Lagrangian  problem  is  sepiirable  into  subproblems  which  are  readily  solved. 

Since  there  are  no  coupling  constraints  between  the  Xpri  and  Zu  variables,  the 
Lagrangian  problem  Dj  can  be  decomposed  into  two  subproblems  as  Zpi{U)  = 
Zd\\{U)  +  Zdi2{U)  where 
Subproblem  Du: 

Zdu{U)  =  min  U,a^/wAV  (8) 

I^L  p^P  reRp 

subject  to 

^  AV  =  1  VpeP 

r€Rp 

0  <  Xpr  <  1  integer 
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and 


Subproblem  Du-' 

Zm2m  =  ^’^T,12iC„-U,Q,)Z„  (9) 

leL  t€T 

subject  to 


Y.Zu  <  I  ^leL 

t€T 

0  <  Zit  <  1  integer 

Subproblem  Du  can  be  separated  into  |  P  \  subproblems,  one  for  each  pair,  where 
the  subproblem  associated  with  the  p‘^  pair  is 


=  mina,  Y1  llZViIr,  AV 

reRp  [iei. 

subject  to 

=  1 

r^Rp 

0  <  Xpr  <  1  integer 

and  then  Z£)u{U)  = 

Subproblem  Du  can  also  be  separated  into  |  L  |  subproblems,  one  for  each  link, 
where  the  subproblem  associated  with  the  link  is 


Zdu{U)  =  mmj;i{Cu-U,Qt)Zu 

i^T 

subject  to 


s  1 

t^T 

^  ^  Zit  <  I  integer 
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and  then  Zoi^iU)  =  JlteT  ^Dui^)- 


E.  IMPLEMENTATION  PROCEDURE 

To  test  the  applicability  of  the  subgradient  optimization  procedure  to  network 
design  problenas  and  to  obtain  estimates  on  the  qucdity  of  bounds  and  solutions 
generated,  this  thesis  implemented  the  procedure  proposed  in  previous  section  in 
turbo  Pascal  on  IBM  PC  computer. 

1.  System  Flow 

The  program  consists  of  two  main  parts.  The  first  is  data  generation,  and  the 
second  is  problem  solving.  Figure  8  shows  the  system  flow.  The  program  is  initiated 
by  entering  the  number  of  nodes.  It  then  creates  data  file  with  source-destination 
node  pairs,  and  their  distance  and  amount  of  data  traffic.  The  user  is  then  prompted 
to  enter  a  number  that  he/she  wantsa  to  generate  a  set  of  cauididate  routes.  Next, 
the  user  is  prompted  to  specify  a  number  of  subgradient  parameters. 

Two  types  of  iterations,  major  and  minor,  are  defined  in  the  program.  A 
major  iteration  consists  of  many  minor  iterations,  each  of  which  is  an  evaluation 
of  the  current  value  of  the  objective  function,  a  choice  of  subgradient  direction  and 
modification  of  the  value  of  U.  The  minor  iterations  terminate  when  one  of  the 
stopping  conditions  is  satisfied  such  as  the  difference  between  the  best  feasible  solution 
and  the  lower  bound  based  on  the  Lagrangian  dual  is  small,  the  subgradient  is  small 
or  an  upper  limit  on  the  number  of  iterations  has  been  reached. 

Throughout  the  optimization  process,  the  program  updates  the  Lagrangian 
value,  the  multipliers  generated  in  each  major  iterations  and  the  best  feasible  solution 
yet  generated. 

Based  on  the  output  generated  after  a  major  iteration,  the  user  may  terminate 
or  reenter  the  subgradient  optimization  procedure  by  changing  the  number  of  routes, 
iteration  counter  or  improvement  counter.  This  iterative  process  continues  until  the 
solution  does  not  get  improved. 


Figure  8:  System  Chart 
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2.  Heuristic 

The  subgradient  optimization  procedure  requires  an  overestimate  Zip  on 
Zd\{U*)  used  to  update  the  Lagrange  multipliers.  This  overestimate  can  be  a  feasible 
solution  generated  by  a  heuristic.  Figure  9  shows  the  heuristic  that  has  been  used  in 
this  thesis  to  generate  good  feasible  solutions. 

In  each  minor  iteration  of  the  subgradient  optimization  procedure,  dual  solu¬ 
tion  and  Z^^  are  obtained.  These  solutions  might  be  infeasible  because  a  set  of 
capacity  constraint  (3)  was  relaxed.  This  constraint,  therefore,  may  be  violated  by 
the  dual  solution  and 

Yi  E^P^riX^>J2QtZit  3/€L. 

p&Pr^Rp  teT 

By  correcting  this  infeasibility  situation  for  all  /  which  has  been  violated, 

X;  E  W€i, 

p€P  T€Rp  t€T 

the  solution  becomes  feasible. 

If  this  feasible  solution  is  smaller  than  previous  best  feasible  solution  which 
is  upper  bound,  it  replaces  the  previous  solution  and  becomes  the  new  best  feasible 
solution. 
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INPUT  :  Infeasible  Dual  Solution  and 
OUTPUT:  Feasible  Solution 

BEGIN 

OBJECTIVE  :=  0; 

FOR  each  Link  I 
BEGIN 
t  :=  0; 

WHILE  >  Qt  DO 

t  :=  t  +  1; 

Zjl  :=  t; 

OBJECTIVE  :=  OBJECTIVE  +  CuZu 
END;  (  FOR  LOOP  ) 

IF  OBJECTIVE  <  UPPER-BOUND  THEN 
UPPER-BOUND  :=  OBJECTIVE; 

END.  (  HEURISTIC  ) 


Figure  9:  Heuristic 
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IV.  COMPUTATIONAL  RESULTS 

The  solutiop  approach  proposed  in  Chapter  III  has  been  tested  using  network 
design  problems. 

A.  Computational  Experiments 

The  procedure  was  tested  under  two  situations:  Under  one  situation  [case  A],  the 
cost  of  T3  circuit  was  9  times  of  cost  of  Tl.  The  other  situation  [case  B]  was  such  that 
the  cost  of  T3  circuit  was  5  times  of  Tl.  Under  each  situation,  we  have  tested  several 
networks  that  have  different  number  of  nodes.  To  find  better  feasible  solutions,  we 
have  tested  this  procedure  several  times  changing  the  number  of  candidate  routes  — 
8,  12,  16  and  20  candidate  routes. 

The  distance  and  amount  of  traffic  for  each  of  source-destination  node  pairs  and 
circuits  information  were  randomly  generated  by  the  random  number  generation  pro¬ 
cedure.  Figure  10  shows  one  of  sample  data  used. 

The  procedure  has  the  number  of  subgradient  parameters.  For  the  test,  we  ini¬ 
tially  set  A®  to  0.5  and  C/®  to  0.  The  stopping  criteria  parameters  are  set  as  follows: 
stepsize  limit  to  0.0001,  error  limit  to  0.005,  and  iteration  limit  to  100. 

B.  Summary  of  Computational  Results 

The  results  of  the  computationad  tests  on  the  five  randomly  generated  problems 
axe  summarized  in  the  Table  1  and  2.  The  values  shown  in  the  tables  represent  the 
feasible  solutions  for  network  with  6,  10,  15,  20  and  25  nodes.  The  asterisk  means 
that  the  previous  feasible  solution  is  not  improved. 

While  testing  the  procedure,  we  had  some  problems  with  the  size  of  main  memory 
and  the  processing  times.  One  of  the  problems  was  memory  overflow  error  when  we 
were  trying  to  enlarge  the  number  of  node,  the  number  of  candidate  route  of  each  pair, 
and/or  the  number  of  circuit.  Since  it  takes  times  to  locate  all  candidate 


30 


node  pair  information 
source  dest.  distance  traffic 


1 

2 

443 

0 

1 

3 

592 

60 

1 

4 

488 

71 

1 

5 

516 

0 

1 

6 

648 

4 

2 

3 

149 

0 

2 

4 

607 

0 

2 

5 

319 

53 

2 

6 

539 

0 

3 

4 

695 

0 

3 

5 

351 

61 

3 

6 

562 

70 

4 

5 

371 

64 

4 

6 

298 

46 

5 

6 

220 

97 

circuits  info, 
capa .  cost 

6 

10 

12 

15 

18 

19 

24 

21 

30 

• 

31 

•  « 

1554 

625 

1560 

627 

1566 

637 

1572 

642 

1578 

646 

Figure  10;  Sample  Input  Data 
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Figure  11:  Graph  of  Time  Complexity 


routes,  we  limited  the  maximum  number  of  hops  of  candidate  routes  to  4  and  the 
maximum  number  of  nodes  to  25.  Figure  11  shows  how  much  time  is  required  as  the 
number  of  nodes  increases. 

Table  1  and  2  summarized  the  results  obtained  for  networks  with  the  different 
Ctindidate  routes.  In  most  of  cases,  the  feasible  solution  improve  when  the  number  of 
candidate  routes  increases. 
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Ni 

Nc 

1 

2 

3 

4 

5 

6 

8 

180175 

253761 

242297 

185416 

1  2 

178883 

» 

238662 

185119 

* 

1  6 

* 

* 

20 

178001 

253413 

* 

183249 

* 

1  0 

8 

629591 

505213 

530912 

613584 

569267 

1  2 

486110 

530421 

* 

1  6 

622310 

449811 

522941 

577069 

568054 

20 

621037 

* 

547983 

1  5 

8 

1415404 

1678792 

1  110691 

1453550 

1373344 

1  2 

1414896 

1 605410 

1060431 

1398207 

1 291821 

1  6 

1388078 

1474272 

1048101 

1382763 

20 

1343018 

•  1 

1355496 

1240806 

20 

2501282 

2645749 

2226695 

1  2 

21620^9 

2227108 

2468560 

2578408 

2148869 

1  6 

2114736 

2176106 

* 

2542664 

20 

* 

• 

2460715 

2509654 

25 

8 

3701524 

3482153 

3497099 

3724747 

4357936 

1  2 

3595955 

3474732 

3302242 

3712604 

4258339 

1  6 

* 

3382492 

3685839 

* 

20 

* 

* 

* 

* 

4176299 

N,  :  Number  of  nodes 
N,  :  Number  of  candidate  routes 
for  each  node  pair 
*  :  Indication  for  non improvement 

Table  1 ;  Summary  of  Fe«isible  Solutions  :  Case  A 
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Ni 

Nc 

1 

3 

4 

5 

6 

8 

127895 

185711 

138250 

171362 

129688 

1  2 

* 

* 

* 

4c 

1  6 

125755 

* 

136375 

4c 

4c 

20 

* 

178148 

134125 

171363 

* 

1  0 

8 

436919 

378817 

378229 

438411 

437988 

1  2 

* 

348660 

* 

417009 

410846 

1  6 

428966 

* 

377508 

* 

20 

326292 

361724 

4c 

4c 

1  5 

8 

1015946 

1184250 

770354 

1000964 

954370 

1  2 

979128 

1156759 

* 

4c 

909834 

1  6 

751964 

973137 

4c 

20 

944654 

1055800 

* 

947506 

894031 

20 

8 

1610162 

1682611 

1859701 

1888310 

1744943 

12 

1500080 

1496639 

1771325 

1749449 

1507031 

1  6 

* 

* 

1682985 

* 

4c 

20 

* 

* 

1708854 

1495353 

25 

8 

2730969 

2678297 

2645483 

2628391 

3166161 

1  2 

2613587 

2578349 

2408956 

3048471 

1  6 

2438277 

2508159 

2391261 

* 

2985384 

20 

* 

* 

* 

2522723 

4c 

Table  2:  Summary  of  Feasible  Solutions  :  Case  B 
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V.  CONCLUSION 


We  have  made  a  mathematical  formulation  for  minimum  cost  design  of  digital  com¬ 
munication  network.  This  thesis  used  the  Lagrangian  relaxation  and  the  subgradient 
optimization  procedure  to  obtain  lower  bounds  and  feasible  solutions  in  minimizing 
the  cost  of  network  design  with  given  node  location  and  traiffic  of  each  node  pairs. 
Computational  experience  indicates  that  this  procedure  is  effective  in  most  instances. 
The  quaJit}’  of  bound  generated  depends  to  a  large  extent  on  the  limits  set  on  the 
improvement  and  iteration  counters.  Low  settings  terminate  the  procedure  before 
it  reaches  a  good  feasible  solution,  while  high  settings  consume  excessive  computing 
resources.  Good  settings  of  these  parameters  require  careful  balancing  between  these 
two  consideration. 

The  difference  between  upper  bound  and  lower  bound  was  fifty  to  seventy  percent 
in  our  experiments.  Our  conjecture  is  that  the  other  constraint  relaxation  is  required 
to  get  a  better  bound.  This  thesis  did  not  explore  that  constraint  relaxation. 

Our  model  can  be  extended  by  considering  reliability  constraints,  which  is  left  for 
future  research. 
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appendix  A:  SOURCE  CODE 


A.  DATA  STRUCTURE  DESCRIPTION 


{*  TITLE  :  Minimuin  Cost  Design  of  Network  *} 

{*  AUTHOR  :  Hwang,  Yong  Goo  *) 

{*  DATE  :  13  June  1990  *} 

{*  REVISE  ;  15  Nov.  1990  *) 

{*  DESCRIPTION  :  This  program  is  to  get  best  *) 

{*  feasible  solution  for  minimum  cost  design  of 
{*  network.  It  has  been  implemented  in  Turbo  *} 

{*  Pascal  and  run  on  IBM  PC  computer.  *} 

{*  System  flow  is  as  follows:  *) 

{*  1.  Generate  data  and  store  it  to  'net.dat'  *) 

{*  file.  *) 

{*  -  node  pairs  information  *) 

{*  -  circuit  information  *) 

{*  2.  Get  data  from  'net.dat'  file  *} 

{*  -  make  route  table  *} 

{*  3.  Solve  the  problem  *) 

{*  4.  Out  report  report 


{******* ***********************************************j 

PROGRAM  NETWORK_DESIGN; 

USES  CRT,  GRAPH; 

CONST 

VLARGE  =  le20; 

MAX_I  =  20;  {  maximum  No.  of  nodes  } 

MAX_P  =  190; 

MAX_R  =  6000; 

MAX_T  =  200;  {  maximum  size  of  circuit  } 

max_nc  =20;  {  maximum  No.  of  candidate  routes  ) 

F_NAME  =  'net.dat'; 

FILE_NAME  =  'out.dat'; 

TYPE 

NODE  =  1..MAX  I; 

PAIR  =  0..MAX~P; 

CNFG  =  0..MAX_T; 

POINT  =  ARRAY [1. .100]  OF  INTEGER; 
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CIRCUTE  =  RECORD 

FT_1  :  INTEGER; 

T_1  ;  INTEGER; 

T_3  :  INTEGER; 

CAPACITY;  INTEGER; 

COST  :  INTEGER; 

END; 

CIRCUTE_TABLE  =  ARRAY[0. .MAX_T]  OF  CIRCUTE; 

SD_PAIR  =  RECORD 

S  :  INTEGER; 

D  :  INTEGER 
END; 

PAIR_TABLE  =  ARRAY [ 1 . .MAX_P]  OF  SD_PAIR; 

VAR 

i;  NODE; 
p;  PAIR; 
t:  CNFG; 

r,ni,np,nr ,nt,nc  ;  INTEGER; 

ORGN , DSTN :  array [ 1 . . MAX_P ]  of  NODE ; 

DISTANCE, TRAFFIC:  array[l. .MAX_P3  of  REAL; 

CAPA, RATE;  array[0. .MAX_T]  of  REAL; 

ROUTE ;  array [ 0 . . MAX_R , 1 . . 4 ]  of  PAIR ; 
px,dx:  array[l. .MAX~P]  of  0. .MAX_R; 
pz,dz:  array [0. .MAX_P]  of  CNFG; 

ITRCTR , ITRLIM , IMPCTR , IMPLIM :  INTEGER ; 

UPPBND , LOWBND , pob j , dob j :  REAL ; 

NORM , ERRORLIM , LAMBDA , LAMBDALIM , STEPS I ZE :  REAL ; 
ALPHA, BEST_ALPHA, SUBGRAD:  array[0..MAX  P]  of  REAL; 
OPTIMAL, TERMINATED, FRAGl:  BOOLEAN; 

DATA_FILE  :  TEXT; 

POINT_X,  POINT_Y  ;  POINT; 

T_TABLE  ;  CIRCUTE_TABLE ; 

S_D_PAIR  ;  PAIR_TABLE; 

ANSI  ;  CHAR; 
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B.  DATA  GENERATION  PROCEDURE 


{*  TITLE  :  Generate  Data  *} 
{*  AUTHOR  :  Hwang,  Yong  Goo  *) 
{*  DATE  ;  15  June  1990  *} 
{*  REVISE  :  15  July  1990  *} 
{*  DESCRIPTION  :  This  procedure  will  make  a  data  *} 
{*  file  that  consists  of  node  pair  informations  *) 
{*  and  circuit  informations  with  random  genera-  *) 
{*  tion  procedure  *) 


I  ****4r***4c*ikikit*********4t*ir*****ift4r4t******it*4r** **********  j 

PROCEDURE  GENERATE_DATA; 

VAR 

t  :  INTEGER; 


I  ******************************************** j 

PROCEDURE  T_TABLE_GENERATION; 

VAR 

i,  j,  k  :  INTEGER; 

BEGIN 

nt  ; =  0 ; 

FOR  i  ;=  0  TO  2  DO 
FOR  j  :=  0  TO  9  DO 
FOR  k  :*  0  TO  3  DO 
BEGIN 

T_TABLE [ nt ] . FT_1  : =  k; 

T_TABLE[nt] .T_l  ;=  j; 

T_TABLE[nt] .T_3  :=  i; 

T_TABLE[nt] .CAPACITY  :=  k*6  +  j*24  +  i*672; 
CASE  k  OF 

1  :  T_TABLE[nt] .COST  :=  10  +  j*21  +  i*219; 

2  :  T  TABLE [nt] .COST  :=  15  +  j*21  +  i*219; 

3  :  T”TABLE[nt] .COST  :=  19  +  j*21  +  i*219; 

ELSE 
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T_TABLE[nt] .COST  :=  j*21  +  i*219; 

END; 

nt  :=  nt  +  i; 

END; 


END; 


FUNCTION  RANDINT(MIN,MAX  :  INTEGER):  INTEGER; 

{*  The  RANDINT  function  suppliers  a  random  integer  within 
the  range  of  the  specified  arguments  min  and  max. 

For  example,  RANDINT (10, 20)  gives  a  random  integer 
between  10  and  20.  *} 

CONST 

MININT  =  -5000; 

VAR 

RANDRANGE  :  INTEGER; 

BEGIN 

IF  (MIN  <  MININT)  THEN 
MIN  ;*  MININT; 

RANDRANGE  MAX  -  MIN  +  1; 

RANDINT  :=  RANDOM (RANDRANGE)  +  MIN; 

END; 


^  4r**4r***4t4:**ir**4t*4r*4t4c4cit4t**  J 

PROCEDURE  NODE_GENERATION; 

VAR 

n  :  INTEGER; 

BEGIN 

FOR  n  :=  1  TO  ni  DO 
BEGIN 

POINT_X[n]  ;*  RANDINT (50, 1000) ; 
POINT_Y[n3  :=  RANDINT (30, 800) ; 
END; 
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END; 


PROCEDURE  PAIR_GENERATION ; 

VAR 

j,  k,  m,  p  :  INTEGER; 

DISTANCES  ;  REAL; 

TRAFFIC_VOLUMN,  TEMP_TRAFFIC_VOLUMN  ;  INTEGER; 
X,  y  :  REAL; 

BEGIN 


p  :*  0; 

FOR  j  :=  1  TO  (ni  -  1)  DO 
BEGIN 

m  ;=  j  +  1; 

FOR  k  :*  m  TO  ni  DO 
BEGIN 


END 


p  ;=  p  +  1; 


S_D_PAIR[p] .S  :=  j; 

S_D  PAIR[p].D  :=  k; 

X  ;=  POINT_X[j]  -  ’■GINT_X[k]; 

Y  :=  POINT_Y[j]  -  POr4T_Y[k]; 

DISTANCES  :=  SQRT((X  *  X)  +  (Y  *  Y) ) ; 
TEMP_TRAFFIC_VOLUMN  :=  RANDINT(-30, 100) ; 

IF  (TEMP_TRAFFIC  VOLUMN  <  0)  THEN 
TRAFFIC_VOLUIW  :=  0 
ELSE 

TRAFFIC_VOLUMN  :=  TEMP_TRAFFIC_VOLUMN; 
WRITELN(DATA_FILE,  j;5,  k;5,  DISTANCES : 10 : 0 , 

TRAFFI C_VOLUMN ; 1 0 ) 


END; 


END; 


{*  main  part  of  data  generation  *} 
BEGIN 

ASSIGN ( DATA_FILE , F_NAME) ; 

REWRITE (DATA_FILE) ; 

RANDOMIZE; 

WRITELN ( '  Enter  number  of  NODES?  ' ) ; 
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READLN(ni) ; 

np  :=  ni*(ni-l)  DIV  2; 

T_TABLE_GENERATION ; 

{*  WRITELN(DATA_FILE,  ni:5,  np:5,  nt-l:5);  *} 
NODE_GENERATION ; 

PAIR_GENERATION ; 

FOR  t  :»  1  TO  (nt-1)  DO 

WRITELN ( DATA_FILE ,  T_TABLE [ t ] . CAPACITY : 6 , 

T_TABLE[t] .COST:  5) ; 

CLOSE (DATA_FILE) ; 

END; 
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C.  SOLVING  LAGRANGIAN  PROBLEM  PROCEDURE 


^^h**************************  **********  *****************  j 


TITLE  :  Generate  source-destination  pair 

AUTHOR  :  Hwang,  Yong  Goo 

DATE  :  18  July  1990 

REVISE  :  18  July  1990 

DESCRIPTION  ;  This  function  provides  a  node  pair 


{* 

{* 

{* 

{* 

{* 

^******************************************************'^ 
FUNCTION  p_(i,j:  NODE):  PAIR; 


*} 

*) 

*) 

*} 

*) 


VAR 

o , d  :  NODE ; 
p  ;  INTEGER; 


BEGIN 


IF  i  <  j  THEN 
BEGIN 
o:=i; 
d;=j; 

END 
ELSE 
BEGIN 
o;*j ; 
d;=i; 

END; 

p:=0; 

REPEAT 

p;=p+l; 

UNTIL  (ORGN[p]  =  o)  and  (DSTN[p]  =  d) ; 
p_:=p; 

END; 


{ ****************************************************** ) 


{* 

TITLE 

:  Get  data  and  Make  route  table 

*) 

{* 

AUTHOR 

:  Hwang,  Yong  Goo 

*) 

{* 

DATE 

;  18  July  1990 

*} 

{* 

REVISE 

:  12  Oct.  1990 

*) 
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{*  DESCRIPTION  :  This  procedure  make  route  table  *) 
{*  from  'net.dat'  file  *} 

{******************************************************} 

PROCEDURE  GET_DATA; 

VAR 

p:  PAIR; 
t:  CNFG; 
datf ;  text ; 

PROCEDURE  BUILD_ROUTE_TABLE ; 

TYPE 

rte  *  RECORD 

11,12,13,14;  PAIR; 
ds:  REAL; 

END; 

ptr  =  -^rte; 

VAR 

i,j,k:  NODE; 

P,lkl,lk2,lk3,lk4:  PAIR; 
c;  INTEGER; 
dists:  REAL; 

ptrs;  ARRAY [ 1. .max_nc]  of  ptr; 

aptr:  ptr; 

nset:  SET  OF  NODE; 

LONGER;  BOOLEAN; 

PROCEDURE  REVISE_PTRS  (k;  INTEGER) ; 

VAR 

c:  INTEGER; 

BEGIN 

IF  k  <  nc  THEN 

IF  dists  <  ptrs[k]^.ds  THEN 
BEGIN 

dispose (ptrs [ nc-1 ] ) ; 

FOR  c;=  nc-1  downto  k+1  do 

ptrs[c] :=ptrs[c-l] ; 
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new(aptr)  ; 
with  aptr'^  do 
BEGIN 

12:=lk2;  13:=lk3;  14:=lk4; 

ds:-dists; 

END; 

pt rs [ k ] ; =apt r ; 

END 

ELSE 

BEGIN 

REVISE_PTRS  (k+1) ; 

END 

ELSE 

LONGER : =true ; 


END; 


BEGIN  {*  main  of  build  route  table  *} 

{*  create  number  of  nc  shortest  routes  *} 

{*  -  with  the  shortest  being  the  direct  link.  *} 

WRITELN('  Wait  a  moment  ii  '); 

FOR  p:=l  to  np  do 
BEGIN 

FOR  c;=l  to  nc-1  do 
BEGIN 

new ( apt r) ; 

aptr^ . ds ; =VLARGE* 10 ; 
ptrs[c3 ;=aptr; 

END; 

nset ; = [ 1 . . ni ] ; 
nset:=nset-[ORGN[p] ] ; 
nset :=nset-[ DSTN [p] ] ; 

FOR  i;=l  to  ni  do 
IF  i  in  nset  THEN 
BEGIN 

lkl;=p_(ORGN[p] ,i) ; 
lk2;=p_(i,DSTN[p]) ; 
lk3 : =0 ; 
lk4 : =0 ; 

dists : =DISTANCE [ Ikl ] +DISTANCE [ lk2 ] ; 
I/5NGER:=false; 

REVISE_PTRS  (1) ; 

IF  LONGER  =  false  THEN 
FOR  j:=l  to  ni  do 
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IF  j  in  nset*-[i]  THEN 
BEGIN 

lk2:*p_(i, j) ; 
lk3:*p_(j,DSTN[p3) ; 

dists : ^DISTANCE [ Ikl ] +DISTANCE [ lk2 ] 

+DISTANCE[lk3] ; 

REVISE_PTRS  (1) ; 

IF  LONGER  *  false  THEN 
FOR  k:*l  to  ni  do 

IF  k  in  (nset-[i])-[j]  THEN 
BEGIN 

lk3:=p_(j,k)  ; 
lk4:=p_(k,DSTN[p3) ; 
dists : ^DISTANCE [ Ikl 3  +DISTANCE [ lk2  3 
+DISTANCE [ lk3  3  +DISTANCE [ lk4  3 
REVISE_PTRS  (1) ; 

END? 

END; 

END; 

ROUTE [ (p-1 ) *nc , 1 3 : =p ; 

ROUTE[ (p-1) *nc,23 :=0; 

ROUTE [ ( p-1 ) *nc , 3  3 ; =0 ; 

ROUTE [  (p-1)  *nc,4  3  ;»=0; 

FOR  c:=l  to  nc-1  do 
BEGIN 

ROUTE[ (p-l)*nc+c,13 ;=ptrs[c3^.11; 

ROUTE [ (p-1) *nc+c, 23 :=ptrs[c3  ^ .12 ; 

ROUTE[ (p-1) *nc+c,33 ;=ptrs[c3^.13; 

ROUTE [  (p-1)  *nc+c,4  3  ;=ptrs[c3"".14; 

END; 

END;  {*  FOR  p  *) 

FOR  c;=l  to  nc-1  do 
BEGIN 

dispose (ptrs [c3 ) ; 
ptrs[c3 ;=nil; 

END; 

END;  (*  end  of  build  route  table  *} 


BEGIN  (*  main  of  get  data  *) 
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ASSIGN (datf,F_NAME) ; 

RESET (datf) ; 

FOR  p:=l  to  np  do 

READLN(datf ,ORGN[p] , DSTN [p] , DISTANCE [p] ,TRAFFIC[p]) ; 

WRITELN('  Enter  the  No.  of  Candidate  routes  for  each  Pair 

READLN(nc) ; 

CAPA[0] ;=0; 

RATE[0] :=0; 

FOR  t:=l  to  nt  do 

READLN(datf ,CAPA[t] ,RATE[t]) ; 

BUILD_ROUTE_TABLE ; 

END; 


{*  TITLE  ;  Find  initial  feasible  solution  *} 
{*  AUTHOR  :  Hwang,  Yong  Goo  *) 
{*  DATE  :  2  Aug,  1990  *} 
{*  REVISE  ;  15  Aug.  1990  *} 
{*  DESCRIPTION  :  This  procedure  gives  us  initial  *} 
{*  feasible  solution  for  problem.  This  *) 
{*  solution  may  be  fully  connected  network.  *) 


^  ic4r4r***il^**** ******************************** *******'****  j 

PROCEDURE  FIND_INITIAL_FEASIBLE_SOLUTION ; 

VAR 

p:  PAIR 7 
t:  CNFG; 


BEGIN 

pob j : =0 ; 

FOR  p:=l  to  np  do 
BEGIN 
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t:=0; 

IF  TRAFFIC [p]  >  0  THEN 
BEGIN 
REPEAT 

UNTIL  TRAFFIC [p]  <=  CAPA[t3 ; 
END; 

px[p] :=(p-l) *nc; 
pz[p] :=t; 

pobj ;=pobj+DISTANCE[p]*RATE[t] ; 
END; 

UPPBND:=pobj ; 

END; 


^ ***************************************************★**} 


{* 

TITLE 

:  Initiate  subgradient  parameters 

*) 

{* 

AUTHOR 

;  Hwang,  Yong  Goo 

*} 

{* 

DATE 

;  13  June  1990 

*) 

{* 

REVISE 

;  15  Nov.  1990 

*) 

{* 

DESCRIPTION  :  At  first  processing,  this 

*) 

{* 

procedure  initiates  sub.  parameters 

*} 

{* 

autometical ly,  and  then  when  user  reenter 

*) 

{* 

the 

system,  he  can  modify  them. 

*) 

I  4r4r4rA4r*4rA*****4t4r4r4c*4r**********ik*x4t*******4r4***********4t  ^ 

PROCEDURE  PREP_FOR_SUBGRADIENT_ITERATION ; 

VAR 

p  ;  PAIR ; 

ANS  :  CHAR; 

BEGIN 

IF  FRAGl  THEN 
BEGIN 

LAMBDA: =  0.5; 

IMPLIM:=  25; 

ITRLIM:=100; 

ERRORLIM: =0 . 0050 ; 

LAMBDALIM: =0.0001; 

ITRCTR:=0; 

IMPCTR:=0; 
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LOWBND:=0; 

FOR  p:*0  to  np  do 

ALPHA [p] ;=0; 

FRAGl  :=  FALSE; 

END 

ELSE 

BEGIN 

WRITELN('Do  you  want  to  change  some  subgradient  parameters?  ') 
R£ADLN(ANS)  ; 

IF  (ANS  *  OR  (ANS  =  'Y')  THEN 

BEGIN 

WRITELN('  LAMBDA  :  LAMBDA  :4:2); 

WRITELN  ( "  Enter  the  LAMBDA?  " )  ; 

READLN (LAMBDA)  ; 

WRITELN ('  IMPROVEMENT  COUNTER  LIMIT  :  IMPLIM  :4); 

WRITELN ( "  Enter  the  IMPROVEMENT  COUNTER  LIMIT?  ' ) ; 

READLN (IMPLIM) ; 

WRITELN ('  LAMBDA  LIMIT  :  LAMBDALIM  :4:6); 

WRITELN ( '  Enter  the  LAMBDA  LIMIT?  '); 

READLN  (LAMBDALIM)  ; 

WRITELN ("  ITERATION  COUNTER  LIMIT  :  ITRLIM  :4); 

WRITELN ('  Enter  the  ITERATION  COUNTER  LIMIT?  '); 

READLN (ITRLIM) ; 

END 

ELSE 

BEGIN 

LAMBDA:*  0.5; 

IMPLIM:*  25; 

ITRLIM: *100; 

LAMBDALIM: =0 . 0001 ; 

END; 

ERRORLIM: =0 . 0050 ; 

ITRCTR:*0; 

IMPCTR:*0; 

LOWBND:*0; 

FOR  p:=0  to  np  do 

ALPHA [p] :=0; 

END; 


END; 


[hie**************************************************** 
{*  TITLE  :  Solve  the  Lagrangian  dual  *) 
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{* 

AUTHOR 

:  Hwang,  Yong  Goo 

*} 

{* 

DATE 

:  21  Sep.  1990 

*} 

{* 

REVISE 

:  11  Nov.  1990 

*) 

{* 

DESCRIPTION  :  This  procedure  computes 

*} 

{* 

a 

Lagrangian  dual. 

*} 

^1t1c1t1e1t1e1e1c1fk’k-k1eic1e1t1cifkifit1fk1fkie1c1t1t1e1c‘k1eit1fk1fk1t1c1e*1i1c1c1cic1c‘k1t1c1c1t1t'^ 

PROCEDURE  SOLVE_LAGRANGIAN_DUAL; 

VAR 

p:  PAIR; 
t:  CNFG; 
r:  INTEGER; 
sum,nin:  REAL; 

BEGIN 

dob j : =0 ; 

FOR  p;*!  to  np  do 
BEGIN 

min:=VLARGE; 

FOR  r;=(p-l)*nc  to  p*nc“l  do 
BEGIN 

sum ; =ALPHA [ ROUTE [ r , 1 ] ] +ALPHA [ ROUTE [ r , 2 ] ] 

+ALPHA [ ROUTE [ r , 3 ] ] +ALPHA [ ROUTE [ r , 4 ] ] ; 

IF  sum  <  min  THEN 
BEGIN 

min : *sum ; 
dx[p] ;=r; 

END; 

END; 

dob j : =dob j  +TRAFFI C [ p ] *min ; 

END; 

FOR  p:=l  to  np  do 
BEGIN 
min : =0 ; 
dz[p] ;=0; 

FOR  t:=l  to  nt  do 
BEGIN 

IF  DISTANCE[p]*RATE[t]-ALPHA[p]*CAPA[t]  <  min  THEN 
BEGIN 

min : ^DISTANCE [ P ] *RATE [ t ] -ALPHA [ p ] *CAPA [ t ] ; 
d2[p] ;=t; 

END; 

END; 
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dob  j  :  =dob  j  +inin ; 
END; 

END; 


{*  TITLE  :  Make  feasible  solution  *} 
{*  AUTHOR  :  Hwang,  Yong  Goo  *} 
{*  DATE  ;  15  Sep.  1990  *} 
{*  REVISE  ;  15  Nov.  1990  *} 
{*  DESCRIPTION  :  This  procedure  will  give  us  *} 
{*  a  feasible  solution  using  infeasible  sol-  *} 
{*  ution.  If  this  solution  is  lower  than  *) 
{*  upper  bound,  then  it  replaces  upper  bound.  *) 


I  1eicicic1c1e1c1fk1i1e‘kie1e1fk1t1i1t"kitit1e‘kic'k1i‘k1e1e1e'k1eic1e1fk1c1eic'k1c1c1fk'k‘k1c'kie'k1e1c‘k  j 

PROCEDURE  PERTURB_FOR_FEASIBLE_SOLUTIONS ; 

VAR 

p:  PAIR; 
t:  CNFG; 
ob j :  REAL; 
flag:  BOOLEAN; 

BEGIN 

obj  :==0; 

FOR  p:*l  to  np  do 
BEGIN 

flag:=false; 
t :  =0 ; 

IF  SUBGRAD [p]+CAPA[dz[p]]  >  0  THEN 
BEGIN 
REPEAT 

IF  t+1  <=  nt  THEN 
t;=t+l 
ELSE 

flag;=true; 

UNTIL  flag  or  (SUBGRAD[p]+CAPA[dz [p] ]  <=  CAPA[t]); 
END; 

IF  flag  THEN 
obj ;=VLARGE 
ELSE 
BEGIN 

dz[p] :*t; 
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Ob j ; =ob j  +DISTANCE [ p ] *RATE [ t ] ; 
END; 

END; 

IF  obj  <  UPPBND  THEN 
BEGIN 

UPPBND :=obj ; 
pz ; =dz ; 
px:*dx; 

END; 


END; 


{*  TITLE  :  Test  subgradient  parameters  *) 
{*  AUTHOR  :  Hwang,  Yong  Goo  *) 
{*  DATE  ;  15  Sep.  1990  *} 
{*  REVISE  :  18  Nov.  1990  *} 
{*  DESCRIPTION  :  This  procedure  checks  that  *} 
{*  terminated  condition  is  satisfied.  *} 


PROCEDURE  CHECK_IMPRO_COUNTER; 
BEGIN 

IF  IMPCTR  >  IMPLIM  THEN 
BEGIN 

LAMBDA ; =LAMBDA/2 . 0 ; 
IMPCTR :=1; 

END; 


END; 

^  4t4t*4c4r*4t****4r***4t*4r****4t4r4t  ^ 

PROCEDURE  TEST_TERMINATION; 

BEGIN 

IF  (LOWBND  >  0)  and  (UPPBND-LOWBND  >  0) 

and  ( (UPPBND-LOWBND) /LOWBND  <  ERRORLIM)  THEN 
TERMINATED : =true 
ELSE 

IF  LAMBDA  <  LAMBDALIM  THEN 
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TERMINATED :  =*true 
ELSE 

IF  ITRCTR  >*  ITRLIM  THEN 
TERMINATED : =true 
ELSE 

TERMINATED : =f alse ; 


END; 


*  h  it  hlcltlt  **11  Itlelt**  Ich  hleleltle  "k -k  *  ********  ********  it  **** -k 


{* 

TITLE 

:  Update  Lagrange  multipliers 

*) 

{* 

AUTHOR 

:  Hwang,  Yong  Goo 

*} 

{* 

DATE 

:  12  Aug.  1990 

*) 

{* 

REVISE 

:  13  Nov.  1990 

*} 

{* 

DESCRIPTION  :  This  procedure  change  Lagrange 

*} 

{* 

multipliers. 

*) 

^******************************************************'^ 

PROCEDURE  UPDATE_LAGRANGE_MULTIPLIERS ; 

VAR 

p:  PAIR; 

BEGIN 

STEPSIZE : =LAMBDA* (UPPBND-dob j ) /NORM; 

FOR  p:=l  to  np  do 
BEGIN 

ALPHA [p] ;=ALPHA[p]+STEPSIZE*SUBGRAD[p] ; 

IF  ALPHA [p]  <  le-4  THEN 

ALPHA[p] ;=0; 

IF  ALPHA [p]  >=  VLARGE  THEN 
BEGIN 

WRITELN( 'WARNING:  MULTIPLIER  TOO  LARGE'); 

WRITELN (dob j : 10 : 2 , NORM: 10:0, STEPSIZE : 10 : 2 , SUBGRAD [p] : 10 : 2 ) 
READLN; 

END; 

END; 


END; 


{*  TITLE  :  Solve  problem  *} 
{*  AXJTHOR  :  Hwang,  Yong  Goo  *} 
{*  DATE  :  7  Oct.  1990  *} 
{*  REVISE  :  16  Nov.  1990  *} 
{*  DESCRIPTION  :  This  procedure  make  a  best  *} 
{*  feasible  solution  using  above  all  subproc-  *} 
{ *  edure .  * ) 


{ ****************************************************** ^ 

PROCEDURE  SOLVE_PROBLEM; 

BEGIN 

REPEAT; 

ITRCTR : =ITRCTR+ 1 ; 

IMPCTR: =IMPCTR+1 ; 

CHECK_IMPRO_COUNTER ; 

SOLVE_LAGRANGIAN_DUAL ; 

WRITELN (ITRCTR: 4, '  Sdobj:10:2,'  ' ,UPPBND: 10: 2 , "  LAMBDA: 10:4) ; 

IF  dobj  >  LOWBND  THEN 
BEGIN 

IMPCTR:«0; 

LOWBND: =dobj ; 

BEST_ALPHA : =ALPHA ; 

END; 

SUBGRAD [0] :=0; 

FOR  p:=l  to  np  do 

SUBGRAD[p] :=CAPA[d2[p] ]*(-!) ; 

FOR  p;=l  to  np  do 
BEGIN 

SUBGRAD [ ROUTE [ dx [ p ] , 1 ] ] : =SUBGRAD [ ROUTE [ dx [ p ] , 1 ] ] +TRAFFIC [p]; 
SUBGRAD[ROUTE[dx[p] ,2] ] : =SUBGRAD [ ROUTE [ dx [ p ] , 2 ] ]+TRAFFIC[p] ; 
SUBGRAD[ROUTE[dx[p] ,3] ] :=SUBGRAD[ROUTE[dx[p] , 3 ] ]+TRAFFIC[p]  ; 
SUBGRAD [ROUTE [dx[p] ,4] ] :=SUBGRAD[ROUTE[dx[p] , 4 ] ]+TRAFFIC[p] ; 
END; 

FOR  p:=l  to  np  do 

IF  abs( SUBGRAD [p])  <  le-4  THEN 

SUBGRAD[p] :=0; 
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NORM:-0; 

FOR  p:=l  to  np  do 

NORM : =NORM+SUBGRAD [ p ] *SUBGRAD [p ] ; 

OPTIMAL;  “true; 

FOR  p;“l  to  np  do 

IF  SUBGRAD [p]  >  0  THEN 

OPTIMAL: “false; 

IF  OPTIMAL  THEN 
BEGIN 

UPPBND;=dobj ; 
px;=dx; 
pz ; =dz ; 

END 

ELSE 

PERTURB__FOR_FEASIBLE_SOLUTIONS  ; 

TEST__TERMINATION  ; 

IF  not  TERMINATED  THEN 

UPDATE_LAGRANGE_MULTIPLIERS ; 

UNTIL  TERMINATED; 

END; 
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D.  REPORT  PROCEDURE 


( ****************************************************** } 


{* 

TITLE 

:  Make  report 

*} 

{* 

AUTHOR 

:  Hwang,  Yong  Goo 

*) 

{* 

DATE 

:  3  Oct.  1990 

*) 

{* 

REVISE 

:  21  Nov.  1990 

*} 

{* 

DESCRIPTION  :  This  procedure  provides  a  report. 

*) 

I  ****************************************************** ^ 

PROCEDURE  DRAW_NETWORK; 

CONST 

MAXLENGTH  =  250; 

MAXHEIGHT  *  150; 

VAR 

CENTER_X, 

CENTER_y  ;  INTEGER; 
i,  j  :  INTEGER; 

TEMP_FT_1,  TEMP_T_1,  TEMP_T_3  :  INTEGER; 
inkey  :  INTEGER;  " 

ckey  :  char; 

DRIVERVAR,  MODEVAR  :  INTEGER; 

^ ************************************** ^ 

PROCEDURE  GRAPHTITLE  (INTITLE  :  STRING) ; 

BEGIN 

SETTEXT JUSTIFY ( CENTLkIEXT ,  TOPTEXT) ; 

SETTEXTSTYLE (TRIPLEXFONT,  HORIZDIR, 4 ) ; 

OUTTEXTXY (GETMAXX  DIV  2,  1,  INTITLE); 

END;  {*  END  Of  graphtitle  *) 


^ **********************] 
PROCEDURE  GRAPHCONTINUE ; 

CONST 

messagel  =  ^CONTINUE?  (y/n) ' ; 
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MESSAGE2  =  'Enter  the  PAIR  number! ' ; 


BEGIN 

SETTEXTJUSTIFY ( CENTERTEXT ,  TOPTEXT) ; 

SETTEXTSTYLE (GOTHICFONT,  HORI2DIR, 3 ) ; 
OUTTEXTXY(GETMAXX  DIV  2,  GETMAXY  -  40,  MESSAGEl) ; 
READLN  (ckey) ; 

IF  (ckey  <>  'N')  AND  (ckey  <>  'n')  THEN 
BEGIN 

SETTEXTJUSTIFY (CENTERTEXT,  TOPTEXT) ; 

SETTEXTSTYLE ( GOTHICFONT ,  HORIZDIR, 3 ) ; 
OUTTEXTXY(GETMAXX  DIV  2,  GETMAXY  -  50,  MESSAGE2) ; 
READLN (inkey) ; 

END; 

CLEARDEVICE ; 

END;  {*  END  of  graphcontinue  *) 


{*******************) 

PROCEDURE  DRAW_ROUTE; 

VAR 

i,j  ;  INTEGER; 

BEGIN 

LINE (10, 10, 10, 450) ; 

LINE(10,450,getmaxx,450) ; 

LINE (cent er_x, 450 , center_x, 10) ; 

LINE ( getmaxx , 10,10,10) ; 

LINE ( getmaxx , getmaxy , getmaxx , 10 ) ; 

FOR  i  ;=  1  TO  ni  DO 

CIRCLE ( (POINT_X[i]  DIV  3) , (POINT_Y[i]  DIV  3),  2) ; 
FOR  j  :=  1  TO  np  DO 
IF  (pz[j]  <>  0)  THEN 
BEGIN 

{*  TEMP_FT_1  :=  T_TABLE [ j 3 . FT_1 ; 

TEMP_T_1  ;=  T_TABLE[j] .T_l; 

TEMP_T_3  ; =  T_TABLE [ j ] . T_3 ;  * } 
LINE(POINT_X[S_D_PAIR[j3 .S]  DIV  3, 

POINT  Y[S_D_PAIR[J) .S]  DIV  3, 
POINT“X[S_D_PAIR[j] .D]  DIV  3, 
POINtIy[S_D_PAIR[J) .D]  DIV  3) ; 
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END; 

LINE(POINT_X[s_d  PAIR[ROUTE [px[ inkey] , 1] ] .s]  div  4  +  center_x, 
POINT_Y [ s_d“PAIR [ ROUTE [ px [ inkey ] , 1 ] ] ♦ s ]  div  4 , 

POINT_X [ s_d_PAIR [ ROUTE [ px [ inkey ],!]]. d]  div  4  +  center_x , 
POINT_y [ s_d_PAIR [ ROUTE [ px [ inkey ] , 1 ] ] . d ]  div  4 ) ; 

IF  ROUTE[px[inkey] ,2]  >  0  THEN 

LINE(POINT_X[s_d_PAIR[ROUTE[px[inkey] ,2] ] .s]  div  4  +  center_x, 
POINT_y [ s_d_PAIR[ ROUTE [px [ inkey] , 2 ] ] . s ]  div  4 , 

POINT_X [ s  d_PAIR [ ROUTE [ px [ inkey ],2]].d]  div  4  +  center_x , 
POINT_Y [ s~d_PAIR [ ROUTE [ px [ inkey ] , 2 ] ] . d ]  div  4 ) ; 

IF  ROUTE[px[inkey] ,3]  >  0  THEN 

LINE ( POINT_X [ s_d_PAIR [ ROUTE [ px [ inkey ],3]].s]  div  4  +  center_x , 
POINT_Y[s_d_PAIR[ROUTE[px[inkey] ,3]] .s]  div  4, 

POINT_X [ s_d_PAIR [ ROUTE [ px [ inkey ],3]]*d]  div  4  +  center_x , 
POINT_Y [ s_d_PAIR [ ROUTE [px [ inkey ] , 3 ] ] . d ]  div  4 ) ; 

IF  ROUTE [ px [ inkey ] , 4 ]  >  0  THEN 

LINE ( POINT_X [ s_d_PAIR [ ROUTE [ px [ inkey ],4]].s]  div  4  +  center_x , 
POINT_Y[s_d_PAIR [ROUTE [px[ inkey] ,4] ] .s]  div  4 , 

POINT_X [s_d_PAIR[ ROUTE [px [inkey] , 4] ] .d]  div  4  +  center_x, 
POINT_Y [ s_d_PAIR [ ROUTE [px [ inkey ] , 4 ] ] . d ]  div  4 ) ; 

END; 


BEGIN  {*  MAIN  PART  *} 
inkey  : =0 ; 

DRIVERVAR  :=  detect; 

INITGRAPH  ( DRIVERVAR ,  MODEVAR ,  "  )  »* 

CENTER_X  ;=  GETMAXX  DIV  2; 

CENTER_Y  :=  GETMAXY  DIV  2; 

GRAPHTITLE ( 'NETWORK  DESIGN' ) ; 

LINE(10, 10, 10,450) ; 

LINE (10, 450, 550, 450)  ; 

LINE(550,450,550,10)  ; 

LINE (550, 10, 10, 10) ; 

FOR  i  ;=  1  TO  ni  DO 

CIRCLE ( (POINT_X[i]  DIV  2) , (POINT_Y[i]  DIV  2),  2) ; 
FOR  j  ;=  1  TO  np  DO 
BEGIN 

IF  (PZ[j]  <>  0)  then 
BEGIN 

{*  TEMP_FT_1  :=  T_TABLE [ j ] . FT_1 ; 

TEMP_T_1  ; =  T_TABLE [ j ] . T_1 ; 

TEMP_T_3  :=  T_TABLE[ j] .T_3;  *) 


LINE(POINT_XtS_D_PAIR[j] .S]  DIV  2, 
POINT_Y[S_D_PAIR[J] .S]  DIV  2, 
POINT_X[S_D_PAIR[j] .D]  DIV  2, 
POINT_Y [ S_D_PAIR [ J  3 • D ]  DIV  2 ) ; 

END; 

END; 

REPEAT 

graphcontinue ; 

IF  (inkey  <>  0)  THEN 
DRAW_ROUTE; 

UNTIL  (ckey  =  'N')  or  (ckey  =  'n'); 
CLOSEGRAPH; 

END; 


^ ************************************** j 

PROCEDXmE  REPORT_BEST_FEASIBLE_SOLUTION ; 
VAR 

p;  PAIR; 
df :  text ; 
i;  INTEGER; 


BEGIN 

ASSIGN (df,FILE_NAME) ; 

REWRITE (df) ; 

WRITELN(df) ; 

WRITELN  (df, 'relative  error  bound:  ’ , 

( (UPPBND-LOWBND) /UPPBND) *100 : 7 : 4 , '  % ' ) ; 
WRITELN  (df, 'LAMBDA;  ',LAMBDA); 

WRITELN (df) ; 

WRITELN(df , 'The  best  feasible  solution:  ', UPPBND; 10: 2) ; 
FOR  p:=l  to  np  do 
BEGIN 

WRITE(df ,p:2, '  [ ' ,ORGN[p] ;2, '-',DSTN[p] :2, ' ]  :  ' 

,ROUTE[px[p] ,1] :4)  ; 

IF  ROUTE[px[p] ,2]  >  0  THEN 

WRITE(df ,ROUTE[px[p] ,2] :4) ; 

IF  ROUTE[px[p] ,3]  >  0  THEN 

WRITE(df ,ROUTE[px[p3 ,3] :4) ; 

IF  ROUTE[px[p3 ,4]  >  0  THEN 

WRITE(df ,ROUTE[px[p3 ,43 ;4) ; 

WRITELN (df) ; 
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END; 
i :  “0  ; 

FOR  p:*l  to  np  do 
BEGIN 

IF  i  <  4  THEN 
BEGIN 

WRITE(df,'['»p:3,']  ',p2[p3:4); 
i : =i+l ; 

END 

ELSE 

BEGIN 

WRITE(df, '[',p:3, ']  ',p2tp]:4); 
WRITELN(df)  ; 
i :  =0  ; 

END; 

END; 

WRITELN(df)  ; 
close (df) ; 

DRAW_NETWORK; 

END; 


^■^r**********************************^ 

{*  main  body  *) 

BEGIN 

FRAGl  TRUE; 

GENERATE_DATA ; 

REPEAT 

GET_DATA; 

FIND_INITIAL_FEASIBLE_SOLUTION ; 
PREP_FOR_SUBGRADIENT_ITERATION ; 
SOLVE_PROBLEM; 

REPORT  BEST  FEASIBLE  SOLUTION; 
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WRITELN('  Try  again?  ;  (y/n)  '); 
READLN(ANSl)  ; 

UNTIL  (ANSI  =  'N')  OR  (ANSI  =  'n'); 


END.  {***  main  PROCEDURE  ***) 

^•k1t1eie1sie1e1e1fkit1c1cie‘k1eifk1t1e1e1eit1c1cit1fk*ie1fk-k1fkifk1c1t1t1tic1e1e'ic1e*ii'^ 
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APPENDIX  B:  SAMPLE  INPUT  DATA  FILE 


A.  PAIR  INFORMATION 


pairs 

sour,  dest,  distance  traffic 


1 

2 

395 

7 

1 

3 

339 

98 

1 

4 

156 

0 

1 

5 

261 

64 

1 

6 

724 

43 

1 

7 

516 

23 

2 

3 

556 

0 

2 

4 

448 

38 

2 

5 

592 

27 

2 

6 

674 

0 

2 

7 

502 

93 

3 

4 

491 

62 

3 

5 

172 

90 

3 

6 

469 

28 

3 

7 

295 

34 

4 

385 

53 

4 

6 

877 

0 

4 

7 

668 

87 

5 

6 

640 

23 

5 

7 

459 

0 

6 

7 

209 

6 
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B.  CIRCUIT  INFORMATION 

CIRCUIT 


capa. 

cost 

6 

10 

12 

15 

18 

19 

24 

21 

30 

31 

36 

36 

42 

40 

48 

42 

54 

52 

60 

57 

66 

61 

72 

63 

78 

73 

84 

78 

90 

82 

96 

84 

102 

94 

108 

99 

114 

103 

120 

105 

126 

115 

132 

120 

138 

124 

144 

126 

150 

136 

156 

141 

162 

145 

168 

147 

174 

157 

180 

162 

186 

166 

192 

168 

198 

178 

204 

183 

210 

187 

216 

189 

222 

199 

228 

204 

234 

208 

672 

219 

678 

229 

684 

234 

62 


690  238 

696  240 

702  250 

708  255 

714  259 

720  261 

726  271 

732  276 

738  280 

744  282 

750  292 

756  297 

762  301 

768  303 

774  313 

780  318 

786  322 

792  324 

798  334 

804  339 

810  343 

816  345 

822  355 

828  360 

834  364 

840  366 

846  376 

852  381 

858  385 

864  387 

870  397 

876  402 

882  406 

888  408 

894  418 

900  423 

906  427 

1344  438 

1350  448 

1356  453 

1362  457 

1368  459 

1374  469 

1380  474 

1386  478 

1392  480 

1398  490 
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1404  495 

1410  499 

1416  501 

1422  511 

1428  516 

1434  520 

1440  522 

1446  532 

1452  537 

1458  541 

1464  543 

1470  553 

1476  558 

1482  562 

1488  564 

1494  574 

1500  579 

1506  583 

1512  585 

1518  595 

1524  600 

1530  604 

1536  606 

1542  616 

1548  621 

1554  625 

1560  627 

1566  637 

1572  642 

1578  646 


APPENDIX  C:  SAMPLE  OUTPUT  DATA  FILE 


relative  error  bound:  63.5411  % 

LAMBDA:  2 . 5000000000E-01 

The  best  feasible  solution:  274700.00 


1 

[ 

1-  2] 

1 

2 

[ 

1-  3] 

4 

13 

3 

[ 

1-  4] 

3 

4 

[ 

1-  5] 

4 

5 

[ 

1-  6] 

4 

13 

14 

6 

[ 

1-  7] 

4 

13 

15 

7 

[ 

2-  3] 

7 

8 

[ 

2-  4] 

1 

3 

9 

[ 

2"  5] 

1 

4 

10 

[ 

2-  6] 

10 

11 

[ 

2-  7] 

11 

12 

[ 

3-  4] 

13 

4 

3 

13 

[ 

3-  5] 

13 

14 

[ 

3-  6] 

14 

15 

[ 

3-  7] 

15 

16 

[ 

4-  5] 

3 

4 

17 

[ 

4-  6] 

17 

18 

[ 

4-  7] 

3 

4 

13 

15 

19 

c 

5-  6] 

13 

14 

20 

[ 

5-  7] 

13 

15 

21 

[ 

6-  7] 

21 

[ 

1] 

12  [ 

2] 

0[ 

3] 

40[ 

4] 

40[  5] 

0 

C 

6] 

0[ 

7] 

0[ 

8] 

0[ 

9] 

0[  10] 

0 

[ 

11] 

16[ 

12] 

0[ 

13] 

40[ 

14] 

16[  15] 

24 

[ 

16] 

0[ 

17] 

0[ 

18] 

0[ 

19] 

0[  20] 

0 

c 

21] 

1 
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